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THEORETICAL INVESTIGATION OF THE RADIATION 
CHARACTERISTICS OF AN ANTENNA. 


By GEORGE W. PIERCE. 


Received May 6, 1916. 


1. Introduction.— For the proper design of a radiotelegraphic 
transmitting station it is important to know the radiation charac- 
teristics of different types of antenna. 

For example, if a flat-top antenna is to be employed, the question 
arises as to what is the best relation of the length of the horizontal 
part to the length of the vertical part, when the excitation is to be 
produced by a given type of generator. It may be known in a general 
way that the greater the vertical length, the greater the radiation 
resistance; it may also be known that the greater the horizontal 
length of the flat-top the greater the capacity of the antenna will be, 
and the greater will be the amount of current that can be made to flow 
from certain types of generator. Now these two quantities, radiation 
resistance and applied current, are both factors in determining the 
out-put from the antenna. 

For a given generator, with known characteristics, the problem of 
getting the greatest output of high frequency energy is a problem in 
the determination of the maximum value of the product of current 
square and radiation resistance of the antenna. 

But this is not the whole problem, for there comes also fnto con- 
sideration the question as to how much of the radiated energy is 
radiated by the horizontal flat-top in what may be a useless direction. 

Again, of the energy radiated from the vertical part of the antenna, 
how much of it contributes to the electric and magnetic forces on the 
horizon, where the receiving station is situated? 

For the solution of these various problems it is important to know 
the radiation characteristics of the antenna in the form of certain 
functional relations. These relations should be known even when 
inductance is added at the base of the antenna for providing coupling 
or for increasing the wavelength to adapt it to the generator. These 
quantities should be known theoretically, since the ordinary measure- 
ments of these quantities do not permit us to distinguish radiation 





ial illest heed eatin inaleatatahs hatin alienate de aleaced ads chee deaee ats we es ne rt - =) 
sade Mata oe tactert Mabe tabatiees Lameiemtadiemen ath adaneal LIE SOLO IS LANG SONNE ONE NE RT YS be Sh IS 


192 PIERCE. 


that is useful from the useless radiation as heat losses and from the 
radiation in useless directions. 

It is the purpose of this paper to give a treatment of this problem. 
Such a treatment is, so far as I know, up to the present entirely lacking, 
but the method here employed is that developed by Abraham ! in a 
very remarkable paper entitled Funkentelegraphie und Elektrodynamik. 
In that paper, Abraham obtained theoretically the characteristics of 
a straight oscillator vibrating with its natural fundamental and 
harmonic frequencies. The present work is an extension of Abra- 
ham’s method to the much more difficult problem of an antenna with 
a flat-top and with added inductance at the base. 


2. Inadequacy of the Conception of an Antenna as a 
Doublet.— Apart from the brilliant investigation by Abraham, all 
other attempts at the treatment of the radiation from an antenna 
assume that the antenna is a Hertzian Doublet. This is only a very 
crude approximation to the facts, for the derivation of the electromagnetic 
field about a doublet assumes that the length of the doublet is negligible in 
comparison with a quantity that is itself negligible in comparison with the 
wavelength. 

Hence, the doublet theory will apply in all of its essentials to an 
antenna, only provided the length of the antenna is not greater than 
one ten thousandth of the wavelength emitted. Of course, it may be 
that at great distances from the oscillator, the theory that it is a 
doublet may not introduce any large errors into certain problems 
such as the propagation over the surface of the earth; but the present 
treatment shows that the doublet theory does introduce large errors 
into computations of such quantities as the electric and magnetic field 
intensities and the radiation resistance of an antenna. It seems 
probable that other problems also should be revised in such a way as 
to replace the conception of the antenna as a doublet by the view of 
it as an oscillator that has a length comparable with one quarter of 


the wavelength. 


3. Method of the Present Investigation.— In the present in- 


vestigation, a doublet of infinitesimal length is assumed at each point of 


the antenna. This is the device used by Abraham. These elementary 
doublets are free from the objection regarding their lengths, as they 
are of infinitesimal lengths, while the wavelength is that due to the 


1 M. Abraham: Physikalische Zeitschrift, 2, 329-334 (1901). 
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whole antenna and therefore is enormously large in comparison with 
the lengths of the elemental doublets. The electric and magnetic 
forces due to each of the doublets is determined at a distance point 
and is summed up for all of the doublets of the antenna, with strict 
regard to the difference of phase due to the different locations of the differ- 
ent doublets. Such a process performed for all points of a distant 
sphere surrounding the antenna gives the total electric and magnetic 
forces at all points on the sphere. Then by integrating Poynting’s 
Vector over the entire sphere, we obtain the total power radiated, 
and from this we compute the radiation resistance and other charac- 
teristics of the antenna. 

The effect due to the vertical portion of the antenna and to the 
horizontal flat-top portion are computed separately, so as to give 
information as to how much energy is radiated with its electric force 
vertical to the horizon and how much parallel to the horizon. 

In deciding as to the proper distribution of the elemental doublets 
along the antenna, the form of the current curve from point to point 
of the antenna is assumed independently. This process is not entirely 
above reproach, because Maxwell’s equations, if they could be properly 
applied to the problem, would themselves give the distribution that 
is consistent with the applied electromotive force at the base of the 
antenna and with the shape and form of the antenna. This step of 
accurately deriving the distribution is, however, at the present time 
not possible of mathematical execution. 

The distribution here assumed for the current in the antenna, as a 
function of the time and of the position along the antenna, is a 
generalization of the distribution assumed by Abraham, and is given 
in the next section. 


4. Assumed Current Distribution.— The form of antenna to 
which the whole discussion is devoted is illustrated in Figure 1, and 
consists of a vertical portion of length a and a horizontal flat-top 
portion of length b. These quantities a and b may have any relative 
values whatever. 

At the base of the antenna is an arbitrary inductance L for varying 
the wavelength. 

The current at any point P’ of the antenna is assumed to be given 
by the equation 


; . Pre . Ir (Xo 
= ~ fes a a= . 
1 = I sin , fsin'y \] ! (1) 
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= length measured along the antenna from the inductance 


PIERCE. 


velocity of light, 

natural wavelength of the antenna without inductance, 

the wavelength with the inductance, 

the current at the point P’, 

to 


the point P’. 
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FIGURE 1. 


The character of the assumed distribution is as follows: The factor 


TC * 7 . « o . 
t means that the current is sinusoidal in time at every point 


d 


of the antenna, with the angular velocity 


2 
sin 


2rc _ 2ne 


or i° (2) 


The meaning of the other factor 


ed 6 (3) 


sin 4 


\ = ! = J (say) 


is illustrated in the diagrams (a), (b) and (c) of Figure 2. 
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If there is no inductance, \ = Xo, and the factor becomes 


) 
J =I cos a (4) 


This is illustrated in (a). 
In the case with added inductance, \ ¥ Xo, and we must keep the 


general form of J given in equation (3). This equation for positive 
values of / gives the upper half of the diagram (b). When/ is supposed 
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FIGURE 2. 


negative the curves obtained continue along the dotted lines of (b) 
and do not give a figure symmetrical with the upper half. To pro- 
duce proper symmetry the absolute value of | must be employed in equation 
(1) when it 18 applied to the distribution of the image to take account of 


reflection. 
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It is also to be carefully noted that when/ = 0, equation (1) becomes 


2 
_ sin —~ t, (5) 


io = I sin DH 7 


so the amplitude at the base of the antenna is 


ss ™ (6) 


Now, finally, when the antenna has a flat-top it is assumed that the 
top part of the antenna is bent over without any significant change in 
the magnitude of the current at the various points. 

When the equation (1) is to be applied to the vertical portion of 
the antenna, we shall call 


[ = 2’, (7) 


where 


z’ = vertical distance from the ground of the point P’ on the 
antenna. 


When the equation is to be applied to the horizontal part of the 
antenna, we shall call 


l=at+r2, (8) 


where 


x’ = distance along the horizontal part of the antenna to any 
point P” on the flat-top. 


The discussion will now be divided into several Parts: Part I. 
Electromagnetic Field Due to Vertical Portion of the Antenna; 
Part II. Field due to Horizontal Portion of the Antenna; Part ITI. 
The Mutual Term in Power Determination. Part IV. Compu- 
tations of Radiation Resistance. Part V. Field Intensities and 
Summary. 
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PART I. 
FIELD DUE TO VERTICAL PORTION OF ANTENNA. 
5. Codrdinates.— Let the origin of codrdinates be at the point 


of connection of the antenna to the ground. Let the z-axis be vertical. 
About this vertical axis as polar diameter, let us construct a system 











FIGURE 3. 


of polar coérdinates in which the position of any point P is given by 
its distance rp from the origin, and the angles @ and @. 


6 = the angle along meridional lines from the pole, 
@ = the angle along parallels of latitude from a vertical plane of 
reference whose position is at present immaterial. 


This system of codrdinates with the positive directions of the angles 
indicated is given in Figure 3. 

If 2’ is the vertical ordinate of any point P’ on the vertical portion 
of the antenna, and r the distance from P’ to P, and if the distance OP 
is large in comparison with 2’, we may write (see Figure 4) 


r= rm — 2 cos 6. (9) 
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6. Field Due to a Doublet at P’.— At a distant point P the 
electric and magnetic intensities due to a doublet of length dz’ and 
charges e and —e at P’ is, by Hertz’s theory, 


a a 7 (t—1/c), (10) 


0 


where 


f(t) = the moment of the doublet 
=e dz’, (11) 
dE, = the electric intensity in electrostatic units, which is en- 
tirely in the direction of 8; that is, of the meridional lines; 
dH@ = the magnetic intensity in electromagnetic units, which is 
entirely in the direction of the parallels of latitude. 
r = distance P’P in centimeters, 
c = velocity of light in centimeters per second. 


The two dots over the f in (10) indicate the second time derivative. 
In writing equation (10), the slight difference in the direction of 
the perpendicular to r from the direction of the perpendicular to 19 








FIGURE 4. 


is neglected in view of the largeness of 7p in comparison with the length 
z’ measured on the antenna. 
Also the r which should occur in the denominator of (10) has been 
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replaced by 7, which can be done without appreciable error for large 
values of r. The same substitution cannot be 
made in the argument of f in (10), for there r 
determines the phase of the oscillation, and this 
phase changes through an angle of 7 for a half 


wavelength, independent of the distance from the 1 -e 
origin. dz 
y +f +-2,€ 


7. Expression of the Field in Terms of 
Current.— We shall next express the moment 
of the doublet and the intensities of the field in 
terms of the current 7 at the point 2’. To do 
this we shall think of the current as delivering a 
charge + e to one end of the element of length 
dz’ and a change — e to the other end of dz’ in a =e 
certain time. A neighboring doublet has a differ- a 
ent current and delivers different charges + e 
and — ¢; partly counteracting the charges of the 
given doublet, and leaving just the charge e — e; that actually occurs 
on the wire. This is represented in Figure 5. 

With this view of the case 





FIGuRE 5. 


and 
7 ® = ede =F ae’ (12) 


Whence, by substituting the value of 7 from equation (1) into equation 
(12) we shall have, in view of (7) and (9) 





whe 9 | 
dE, = dH, = == Se = (ct — my — 2’ c088): 


in = (% o\ al 
sin \~ M2) de (13) 


By integrating this expression from z’ = 0 to z’ = a, we obtain the 
electric and magnetic intensities at the point P due to direct trans- 
mission from the vertical portion of the antenna. Indicating this 
integration, we have 











200 PIERCE. 


2a I sin 8 
ACI 


a ) . 
Eo = Hg = { cos % (ct — ro — 2 cos 9) 


2r Ao , ’ 
sin) ¢ ‘| dz’. (14) 


By reflection from the earth, which we shall regard as a perfect 
reflector, we have intensities that must be added to the above. These 
intensities may be obtained by considering the radiation to come from 


an image point at a distance z’ below the surface. The effect of this 


is obtained by changing the sign of the z’ in the cosine term of equation 
(14), but as was pointed out in section 4 the sign of z’ in the sine term 
must remain. We obtain thus for the intensities due to the reflected 
wave emitted by the vertical portion of the antenna the value 


_ 2x1 siné 


"a 2r ; 
ees J cos +- (ct — 7 + 2’ cos @)> 


in? (2) ae 5) 
sin = (5 z) dz’. (15) 


Adding the equation (15) for the reflected intensities to the direct 
intensities of (14), remembering that if A and B are any two angles 


cos (A — B) + cos (A + B) = 2 cos A cos A, (16) 


we obtain for the total intensities at P the equation 
: 4r I sin 0 2r " Qa’ 
Eo = Hg = cos — (ct — 1) cos cos 6 }« 
Acro e/a 


. 2 Xo \ / — 
sin = — 2} dz, (17) 


which resolves into 


I sin 6 2 
Ee= H,= or i om COS * (ct — ro) { 


Acro 
. Io {’ 272’ cos 6 2rz Ie! 
sin —- cos cos ——- dz 
2 e/7o PN 
a 2rz’cos6 . 2rz ,, 
— cos cos sin dz (18) 
2 e/a IN r 
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This expression may be integrated by the formulas 360 and 361 
of B. OwPierce’s Short Table of Integrals and gives 





edited = _ ie = (ct — ro) 
cos B cos (A cos 6) — sin B cos 6 sin (A cos 6) — cos G (19) 
where 
9 
D = aria 
9 ; 
| A= “ | (20) 
| G= ~ =A+B | 


The quantity b, which is the length of the flat top, gets into (20) 
' and (19) by reason of the fact that a + b = the whole length of the 
antenna, so that 


No = 4 (at BD). (21) 
Z 
P 
Yo 
' 
= 
a\¢é 
O 





FIGURE 6. 


Equation (19), with the notation of equations (20) and (21) is the 
gencral equation for the electric and magnetic Intensities at any point P, 
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due to the whole vertical part of the antenna. In this formula, referring 
to Figure 6, ° 


ro = the distance OP in cm., 

6 = the zenith angle ZOP, 

b = length of the horizontal flat top in meters, 
length of vertical part of antenna, in meters, 


a= 

Ay = 4 (a + Db) = natural wave length in meters, 

\ = wave length in meters actually emitted, and differing from \o 
by virtue of the added inductance, 

I, = amplitude of current in absolute electrostatic units at the base 
of the antenna and related to I by the equation, 

[y= Tsin =. 


We shall reserve comment on this equation until after investigation 
of other characteristics of the radiation. See Part IV. 


8. Total Power Radiated from the Vertical Part of the 
Antenna. — Having obtained in equation (19) the electric and mag- 
netic intensities at any required point at a distance from the antenna, 
we shall next compute the total power radiated from the vertical part 
of the antenna, and shall then obtain its radiation resistance. 

Since £, and Hg are perpendicular to one another and perpendicu- 
lar to mo, we have, according to Poynting’s theorem for the power 
radiated in the direction of ro through an element of surface dS per- 
pendicular to ro the quantity 


c 
| —_ yi S. 22 
dp 4, bos ds ( 


Let the element of surface be an elemental zone on the surface of 
the sphere, then 


dS = 2nxr,* sin 6dé (23) 


This quantity, together with the values of Kg and Hg from (19), 
substituted in (22) and properly integrated, gives for the total power 
radiated through the whole hemisphere above the earth’s surface, 
the value in ergs per second following: 











f 








Pex 6s . 
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Q]J2 2, a/2 32 (4 . 
p= ai cos” 3 (ct — 10) | cos B f = — a 











( sin 6 
eaien f ™/2 eos? 6 mt - cos 6) dé —_— af" 2 
— 2sin Bcos B cos @ sin (4 = ’ cos (A cos 6) dé 
— 2cos BeosG | */? cos 0s 6) dé 
+ 2sin B cos G f yy cos @ ain (ons 6) “|: (24) 


This equation when integrated gives the power radiated from the 
vertical part of the antenna. The integration is a tedious operation, 
and is given in the next section, which may be omitted by readers 
not interested in the mathematical processes involved. The result 
of the integration is found in Section 10. 


9. The Integration of Equation (24). — By the use of such trigo- 
nometric equations as 
. 1 + cos 22 
cos? x = = 


— 1 — cos 22x 
sin* z = - ; 


the squares of sines and cosines in the integrands of (24) may be 
avoided, and equation (24) written 


or ) an ) 1 iets : dé 
p= yy — + (ct — ro) \ (5 + “os a) ff sin y] 


x/2 : SI 2 a/2 . 
4 08 2 or cos (2 A cos 8) dO at ad B { sin 6 d@ 








sin @ 
sin? 2 B 
— sin 6 cos (2 A cos 6) dé 
~ vo 
sin 2 B & cos 6 sin (2 A cos 6) dé 
2 sin 8 
. _ {77 cos (A cos 6) dé 
— 2eosBeosG 
Jo sin 6 


ax/2 
+ 2sin B cos af ome sin ( (A — a) “|. (25) 


sin 0 
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The third and fourth terms may be integrated directly. In the 
other terms let us introduce a change of variable as follows: 


Let 


u = cos 6 
— du 

dé = ———» 
9 sin 6 


then 


®*7/2 dé *y eons du | [ ] 1 
* = ————4 - ws / 
J sin 6 J lL~ 2 e O (; - u T l—u ne 


1 (" du =k du | {> du vias 
=o, ite’ 3J.114¢s "8s. i+ © 


With this operation as a model, two of the other integrals of (25) 
may be written, respectively 





{~ COs (2 Acos6)d6_ 1 f° COs (2 Au) du (27) 
nS sin 0 2 J-1 i+e fe 
. cos (A cos6) dé I i cos (Au) du (28) 
Jo sin 6 Se Il+u - 


Another of the integrals, examined in more detail, gives 


r cos 6 sin (2 A cos 6) dé 


sin 6 


‘fk usin (2 Au) du 
vo 


l—w 


: {~ ( a. sin (2 Au) du 
em ing ital” 
>. { sin(2Au)du . 1 eg sin (2 Au) du 





1+ u = We 1+ u 
ae fo sin (2 Au) du . 
Baus It+u * (29) 


Similarly, the remaining integral becomes 


sin 6 o a 


[~" cos@sin (A cos6)d6—— sin (Au) du 
e/J0 =~ I 7 u 








ee Er, OE Cw ae 
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Returning now to equation (25), we shall integrate the third and 
fourth terms, setting them first, and shall substitute (26) to (30) for 
the other terms, obtaining 





212 .\ Qe / sin? B . sin? Bsin2 A 
P = 3 COs” > (ct =< ‘ | - -— oe ac 94 
1 > wr 1 ig du 
+ (5 tos a), -. 1l+u 
4. 008 2 2B 22 +! cos (2 Au) du 4 sin2B (‘*! sin (2 du) du 
1+ u So as 1+ u 
' tT! cos (Au) du ‘+1 sin (Au) du | 
—_ S - " d ) 
cos G bail ‘ee + sin By ey | (31) 


Let us now write 


y¥y=2A(14+ wu) 
2 Au=y — 24, 
dy 
du = a4 
du dy. 
1+ u an 


then the second and third integrals of (31) become 


cos2B ('*! cos (2 Au) du , sin2B [('*' sin (2 Au) du 
4 . 1+ u ) 4 J —1 1+ 


cos 2B ly 
= - - } {cos y cos 2 A + sin y sin 2 A} 
a eJo ° 








sin 2 B o., , d+ 
a | {sin y cos 2 A — cos y sin 2 A} - : 
4 Jo x 
_ cos (2A + 2B) [~ cos y dy 4 sin (2.4 + 2 B) 'g sin y dy 
4 eo 7 4 Jo Y 


cos2G (‘44 cos sin2G [Asin y 
6 See | 2 dy + - | : dy . 
4 eo Y + e/o 7 





In like manner, the last line of (31) becomes 


A °2A a’ 
— COS in ¥ ita 
— cos’?G { ; a dy — cos G sin G | mae dy. (33) 


Ay 
I 
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Let us now decompose the coefficient of the first integral of (31) as 


follows: 
1 , cos?G 1 ,  ¢os? G 
/Ai tl wie 
1 ] s2G 
“7° + ‘7 + cos’? G 
= _ cos2 Se 4 cos? G. 


Then the whole equation (31) may be written 

















| 20°. 7 om Oot ra) sin?B , sin? Bsin2 A 

= - Os* - i -|— —+ — 

a ) d 7 | 2 4A 

cos2G (44(1—cosy)dy , sin2G (‘*4siny 
— = dy 
0 7 
‘ °2A (1 — cosy) d sin2G (?4sin : 
+ cos? G | 2a | al 
Oo 7 - Oo Y 


The various integrals may now be obtained by expanding in series 
and integrating term by term. This gives 


oT } 2 9 
p =~ “cos? \ 20 (et 1») | the aCe A - 1) 





2 2A 

Z dias _ Ay, AP _ | 

4 ) 212 4!4 6!6 

4 L+cos 2G { (24)? _ (2A, 2A _ | | 

2 ) 212 414 6!6 

sin 2G{ ,, _ (44, 44)e_ 

TG] } - 313 T 515 
_sin2G{ ,, _ (24) , (24% : 
at Bee 313 515 | (35) 


Let us now eliminate B from the first terms of this equation, by 
substituting B = G — A, obtaining 











Re ee dias 
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sin 7B /(sin 2A _ 1 — cos 2B /sin 2A ' 
| 94 1) = 2. 3. 
(1 cos (2G — 2A) } /sin2A 
oS ae ‘ — 
| 4 4 \\ 2A 
_s cos 2G cos 2A sin 2G sin 2A 
4 4 4 
sin 2A __ COS 2G sin 4A 
SA 16A 
sin 2G 1 — cos 4A 
hn ion Rta uxcoel 3 
4 4A (36) 








If now we expand in series the quantities involving A in (36) and sub- 
stitute in (35), we obtain, if 











k = 2A 
37 
q = 2G t ) 
_ 22 ge) 27 ({1$ _ ew we 
lie taal ae nd ta 31° BI 71 
Ls 2k 2k! okt 
4) 212° 414° 616 
cos q J a, 
7a t* “Sa” ar 4 
sin ¢ | k tee 
77 1" “arta Rt 
_ cos q J 2k? (2k | 
-—- ae 
_ sing) 2k _ (2k) | (2k) _ om 
4 ) 2! 4! 6! | 
_ cosg§ (2k)? (2k | (2k) _ 
4 ) 212 414 " 616 
csq{ @ eM, el 
T > ) 29 ~ 447 68 
sing§ 5 _ (2k) , (2k) _ | 
re ” 313 515 
sin q { ie k° 
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If now we add together the terms multiplied by sin g and those 
multiplied by cos g, and those not involving g, we have (on factoring 
out the +) 








a : ail ) : ee ik ‘ ote = ed cry r | 
i a eee 
6 —s 8 ys 4 8 - 10 aa | 
+ sing, — 3 a —- 3 + a ta 1s 
Z+2—-9, g4+2—11.. ) : 
| 8!7 aie a SHB mo : 


(39) 


Equation (89) gives the total power radiated by the vertical portion of 
the antenna into the hemisphere above the earth’s surface. In this equa- 
tion, the current factor J is in absolute ec. g. s. electrostatic units, and 
the power p is in ergs per second. 

It is convenient to change the current factor into amperes and the 
radiated power into watts, which can be done by multiplying the right 
hand side of (39) by 30c. This is done, and the equation is rewritten 
in the next section 


10. Result of the Integration for Power Radiated from the 
Vertical Part of the Antenna. — By equation (39), when reduced to 
practical units, the total power radiated into the aérial hemisphere 
from the vertical part of the antenna may be written 


— oe ae 
p = I cos’ > (ct — 7) ( hen — R,cos gq — R3 sing |, (40) 


where 
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1) 2t2 442, , +2, 
fy = 19) “315 . sa” + 76 * | 
a | 
\ 
2 2 2 — ... £2 
R, = 15} = 5 : k? — et 
e+ 25 — 8 PY bn 
7! 2 
~(#+2—5 F+2—7,, | 
R3 = 15} 413 i — 615 k + 
oe. | 
+ i 9 Lt — i | 
Sid \ ) 
q= | 
. . (42) 
= dra 
r 





a = length of vertical part of antenna in same unit as X 
(e. g. meters), 
p’= radiated power in watts instantaneous value, 


fa ——., (43) 


where 


Io = amplitude of current at the base of antenna in amperes. 


11. Radiation Resistance of Vertical Part of the Antenna. 
— In equation (40) is given the power radiated from the vertical part 
of the antenna, on the assumption that radiation from the horizontal 
part of the antenna does not interfere with it. It will be shown later 
in §14 et seq. how this interference is computed and allowed for. 
Accepting for the present the assumption of non-interference, we may 
obtain the radiation resistance of the vertical part of the antenna. 

The radiation resistance is defined as the time average of radiated 
power divided by the time average of the square of the current at the base 
of the antenna. 

In taking the time average of the power (40), it is to be noted that 
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\ 2r ia ae 
the time average of cos? - (ct — 70) cis 1/2. The time average of 


re oe 
current square at the base of the antenna, by (1) is 9 I? sin’ ~~ 


— q — ; 
I? sin? (?) Whence the radiation resistance becomes in ohms 


= R, — Re cos ;- Rs; sin q ’ (44) 
sin’ (2) 


in which R,, Re, R3 and qg have the values given in (41) and (42). 
We shall later give tables of R,, Re, and R3, that will reduce the 


calculations of R to very simple operations, and shall compare the 
results with calculations on the doublet hypothesis and with observa- 


Rg = 


tions. 

We shall, however, first investigate theoretically the radiation from 
the horizontal part of the antenna. This is a problem of considerable 
mathematical difficulty but is capable of solution. 
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PART II. 


FIELD DUE TO HORIZONTAL PoRTION OF ANTENNA. 
12. Introductory Notions.— To determine the electromag- 
netic field and radiation characteristics of the horizontal flat-top 
portion of the antenna, let the rectangular coérdinates of any distant 
point P (Fig. 7) be 2, y, z 


Pp 














P 


FIGURE 7. 





** e , ° , . 
And let the codrdinates of any point P’ on the flat-top be 2’, 0, a; 
the codrdinates of the image point P” be x’, 0,—a. 


Then the distance from the origin of codrdinates to the distant 
point Is 


OP 





mn=Ve+y+? 


The distances of the distant point from the point on the flat-top 
and its image respectively are 
PP=r=Vx—o P+ yt (2 — a)’. 


and PY” P= rr" = V(x 


—xr)?+ y+ (2+ a). 
f 9 9 \9 | 9 9 2 
Then r — ro = Vr — a’ )P + P+ (ec — a)? — Ve + y+ mS. 
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As an approximation, let us multiply by the sum of these radicals and 
divide by the approximate value of this sum for large values of 1; 
namely, by 270, obtaining 


yg’? — Qax’ — 22a + a? 








r = m + —— ae (45) 
fs 4 / 9 2 
r’ =m+ & _. 2za + a ; (46) 
2ro 


13. Determination of Electric and Magnetic Intensities due 
to Flat-top. — The values of r’ and r” in (45) and (46) may be re- 
placed by 7 in intensity factors, but not in phase terms, and give 
for the sum of the effects of a doublet at P’ and another at P” (the 
image doublet) the electric and magnetic intensities 


i, ible ee mt rbe~v'0 , (47) 


roc? ( 


where fi(t) and fo(t) are the moments of the two doublets respectively. 








ne 
P/ JAZ 
Y 
y 
f ic! 
° » 
O x, ¥ 
FIGURE 8. 


The angles ¢g and = correspond to the angles 6 and ¢ of Figure 3, except 
that the figure is turned on its side, so as to put the polar diameter 
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along the z-axis instead of the z-axis. This arrangement is shown in 
Figure 8. The plane of the zero value of = is now to be fixed as the 
plane of the x and z-axes. 

Now using the current distribution of equation (1), we must replace 
l by a+ 2’, which gives, when treated as (12) was treated, 





:;_ 2ael { 2m ( vais x’? — 2z’x — 22a + -) ( 
1 X Cos r C ro Or 
. 27 Ao en fl ( ’ 
sin) ( Credited i dx’. (48) 


The fictitive current at P” is just equal and opposite to that at P’, 
with, however, a different distance from the point P, so we may write 


<= 2enl + —— x? — 2ex! + 22a 4 =) ( 
rs rN wall @ . 2ro \ 
* \ 2m (Xo ’) (7. 
sin ht a a— x) dz’. (49). 


Whence by addition, employing the trigonometric relation 
cos (a + B) — cos (a — B) = — 2sinasin8, 


equation (47) becomes 


4nIT sin . \ 27 xv’? — 2re’ + a*\ 
= dHs = — sin - t—rn— - ) 
dE, = dHy oak sin : = (ct r. - }i 


sin anes in ) x(F } la ‘ 
ed BY Dats 2 al | 


Pra 
— by J 


In this equation we may as usual replace - x \. Also we may 


make an approximation as follows: For large values of ro 


x? — 22a’ + a? rr’ ' 
a =— —-=-2 cosy. 
<1 "0 
aD) » 
' , Hoa ined . 7+ 2 | 
In making this approximation the neglected term is —>—— , and 
<To 
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this is to be neglected even in the phase angle, because its value is 
absolutely small. We have then 


: 4rI sin , @z, 27 
dk, = dHy = — aah v sin 4 sin ) \ (ct — ro + 2’ cos ) 


. | 2rfNr | oe nie 
‘sin ) =(3 —_-a_- x’) ( ix’ | (50) 


This equation may be shortened up by writing 


9 
r= = (ct — 10) (51) 
and 
2a No oarh 
—_— sae —_ — = 52 
B x (> a) x ( 


To obtain the total electric and magnetic intensities due to the 
flat-top, the equation (50) must be integrated for all the doublets 
and their images between the limits 


x’ = 0 and 2’ = 0 


where is the length of the flat-top. This integration is expressed 
in the following equation. 


; 4nIsiny . a fs / 2rrx’ P 
Ey = Hy = — rch sin a sin d a \ COs v) sin 
: esal 
(B - — ) ar’ (53) 


To perform the integration let us introduce a change of variable 
by putting 


Iara’ r 
s=B- \ then dz = — a ds 


and the limits of integration become 


forz’ = 0, s= B, forz’ = 6, =: = 0. 


/ 


Equation (53) then becomes 








yee 











SRO 


ros 





1 2 eh 
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' 2Isiny . Az [(°. , 
E, = Hs = —— sin — sin (r + Beosy — scosy) sins ds 
ro B 


roc 


2I[siny . Azl. . ; 
= ———- sin —| sin(r + Beosy) cos (s cosy) sins ds 
roc ToL JB 
*o 
— cos (tr + Beosy) | sin (scosy) sin sds }- (54) 
JB 


The expressions of this equation may be integrated by the use of 
formulas 360 and 359 of B. O. Peirce’s Tables and give 


2i060Cli«wsaéaHa | ) : 
—— sin - sin (r + Becosy)| — cos s cos (s cos y) 

roc sin ro | 

) 


— cosy sin s sin (s cos ¥) | 


B 


— cos (7 + B cos y) cor y sin 8 cos (s cos y) 


{) 
— cos s sin (s cos Y) | 
B 


2i 


. Az] . 

—_—§— sIn -| Sin (r oa B COS y) ) — | + COS B COS 
roc sin Y ro ( 

B (cos ¥)) 

; ; \ 

+ cosy sin B sin (B cos y) ( 

+ cos (r + B cosy) cos ¥ sin B cos (B cos ) 


— cos B sin (B cos y) | 


| ia \ | 
= . sin — | sinz + cos B — cos (Beosy) - 
roc sin Yo 


+ cosT cos Y sin B — sin (B cos py) i (55) 


Equation (55) gives the electric and magnetic intensities due to the 
flat-top at any distant point whose codrdinates are 


ro = distance of the point from the origin, 

z = vertical height of the point above the earth’s surface, 

y = angle between 7 and the x-axis; this x-axis being parallel to the 
flat-top. 
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The quantities A, B, and 7 are defines by equations (20) and (51). 
We shall next discuss the total power radiated from the antenna. 


14. Concerning Power Radiated from the Total Antenna.— 
It is to be noticed that the electric and magnetic intensities due 
to the flat-top of the antenna and those intensities due to the vertical 
portions of the antenna are directed along the meridional and lati- 
tudinal lines of two systems of polar coérdinates with their poles one 





FIGURE 9. 


quadrant apart. This does not make the respective intensities 
perpendicular to each other, and it becomes necessary to resolve one 
set of these intensities along and perpendicular to the other set of 
intensities. At a given point on the sphere about the origin of co- 
ordinates, the quantities ¢, 6, 2, and y are oriented in a manner 
represented in Figure 9. 

If we let 

a = angle between y and 6 

then also 


angle between @ and 2%. 


~ 
| 


It is also apparent that 


T 


) 


Angle between 2 and @ = a — 


OT 
Angle between y and ¢ = > — 4 


Let us now resolve Ey and Hs into components along @ and per- 
pendicular thereto (that is, along ¢) obtaining for the 6-components 








ais BT da alae 5 ‘ tee 7 7 aa Re ie penta ole OF sera aE $ P 
cane aT ee eae ode och LTA eta Mat 




















Pe ae ca as eel ree Sik. Rae 
© * a a. LY Ae daar 














RADIATION CHARACTERISTICS OF AN ANTENNA. 217 
Ey,.0 = Ey cosa 


= OS —_ ™ \ = > , 
Hys.¢ = Hs cos (« 2) Hs sin a. 


and for the ¢-components 


9 


-_ 


; , oT — 
Eye = Ey cos ( —a})= — £ysina 


Hys,, = Hy cos a. 


Adding these quantities to the corresponding components of the 
intensities due to the vertical part of the antenna, we obtain for the 
total intensities, which are designated by primes, the values 


E", — | OF + ky COS a, 


E’, = — Ey sina, 
H'’,= #£Hysina 
H’, = H, + Hz cosa. 


All of these intensities are perpendicular to 7p. To get the power 
radiated through an element of surface dS perpendicular to 7, we may 
make use of Poynting’s vector, in the form 


dp = 3 ( : H’) dS 
4ir\ / 


where the cross between the vectors means the vector-product. This 
vector-product, expanded, gives 


dp a : (E's _- saa E'. H's) dS 
4dr 


” ra (Fi He + Ey Hs cos” a —- H, ky cosa a Es Hs cos a + 


Ey Hy sin «) ds 


= = (F He + Ey Hz + 2 cos a Ey H, } dS. (56) 


We have already found the first term of this power and have ob- 
tained its integral all over the aérial hemisphere. This integral we 
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have called the power radiated from the vertical part of the antenna. 
We shall call the second term above (56), when properly integrated, 
the power radiated from the flat-top. The third term, since it contains 
both sets of codrdinates, may be called power radiated mutually. 
These designations are merely for convenience in paragraphing the 
mathematics involved. 


15. Power Radiated from the Flat-top.— Let us now enter 
upon a determination of the power contributed by the second term of 
the right hand side of equation (56), and integrate this term over the 
aérial hemisphere; that is, the hemisphere above the surface of the 
earth regarded as a plane. 

The element of area of this hemisphere is 


dS = re sin y dy dd. (57) 


This is to be substituted in the required term involving Ey and H3; 
but these quantities involve the coérdinate z, which must be replaced 
by its value in polar coérdinates 


2 = 7 siny cos 2. (58) 


Besides (57) and (58) we are also to substitute the values of Ey and 
Hs from (55) into the term 


E, and Hy are identical, by (55); the product will give certain 
terms involving sin’7, other terms involving cos’7, and still other terms 
involving sin 7 cos 7; where 7 has the value given in (51). If we take 
the time average for a complete cycle, or, if we prefer, for a time that 
is large in comparison with a complete period, we have 


av. sin’7 = av. cos’r = 5; 
while the average of the product 


av. sin7t cos Tr = 0. 


The integral form of (59) then becomes, if p = the time average 
of radiated power, 
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nde oh A 2d aint — | | 
p= | al cos* B + cos*y sin? B + 1 — 2 cos B cos (Beos y) 


r/2 


—"" | i ‘ 
— 2cosy sin B sin (B cos y) | dX |sin? (A sin y cos D)| | 
. —7/2 ; 





(60) 
We shall first perform the integration with respect to 
4/2 ~ | 
, = | 1 — cos (2A siny cos> 
fa {sin (4 sin cosz)} = fe ‘ s ( v : 
—x/2 a = 2 
2 
™ 
r 1 [(* ’ , - 
a. ms os (2A siny cos 2) dd. 
3 
x /2 
er 1 f° ; 1 - 
= 5 —5 | cos (24 siny cosz) dz — 5 | cos (2 A siny cos >)d> 
- - — - 0 
(61) 
S22". « , - . 
at toot cos (2 A sin y cos 2) db. (62) 
_ = @ Oo 


This last step consists in changing the variable of the first integral 
of the right-hand side of (61) by putting 


w/ feo ww 


which makes the limits i and z+ without any other change, except 
the change of = to >’. But since this is the variable of integration, 
the prime may be omitted, and the terms of (61) added, giving (62). 

Equation (62) may now be integrated by Formula (11) Art. 121 
of Byerly’s Fourier’s Series and Spherical Harmonics giving for the 
integral of (62) 
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~ 
| d= \sin? (A siny cos=)} = —<Jo(2Asiny), — (63) 


to 


where Jo is the Bessel’s Function of the zeroth order, with a develop- 
ment of the form 
ro x og 


iG) = 1-=+—.+ 


2 24 2446 





7 (64) 


Before substituting in (60) let us simplify the general trigonometric 
factor in the brace of (60) by placing cos*¥ by 1 — sin’, and letting 
i; = 2A, as in (42), we then obtain 


2 — sin*y sin? B 


a [° { \ 1 — Jo (k sin y) ( | 
l 4eJo | sin \ | 


— 2cosBcos(B cos y) — 2 cosy sin B sin (B cos ) dy 


_P (*\Rsin?y Fk sin’y , ksin’y 


~ 4edo ) 2 24 2? 4? 6? 
2 — sin’ y sin? B — 2 cos B cos (B cos y) 
elias 
— 2cosy sin B sin (B cos y) ( dy. (65) 


or 


: IP - 2s }.” Fr il 
p= 4.| —22(— 1? spenae J sin’ edd 


i 
| n i” or 

+ sin? BY (— 1)? acronis | sin™+! Ydy 

2°4°0*: > n° eJ0 

n ie” "x . 
+2cos BY (— 1)? ae | sin” y cos (B cosy) dy 

ds ee“ ae 
}.” T 
= { sin”! Y cos W sin 

MW” eo 


(B cos y) dy ] 
n=2, 4, 6,..... (66) 


+ 2sin BY (— 1)2 


224262... 
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Treating these several integrals separately, we have 
T 
7 eS or 
f sin"! Ydy =| sin”?! yYdy +f sin”? Ydy 
Oo e7O T 


a 
ri 


[sive Ydp + [ ‘cos Ydy 


>. 


\ 2-4-6---:n —2/ 
— ). - 
~ 4} 1°3-5 n—l1 } (67) 
by B. O. Peirce’s Tables, Formula No. 483. 
Likewise 
e \2:4-6-+--n 
sin"*+! Ydp = 2- = 
J wer (1-3-d°-ssn +1) (68) 


Now by Byerly’s Fourier’s Series and Spherical Harmonics Equation 


(9), Art. 122, 


| sin” W cos (B cos yy) dy = i J n-i (B) (69) 
vo B a & 


where J n-; (B) isa Bessel’s Function of the order (n — 1) /2, and 


yy 


_ 


jf . ‘ , an 
T (5) is the Gamma Function of = 


For the last integral of (66), we have by Problem 2 and Equation 
(9) of the same article of Byerly’s Fourier’s Series 


"x 
| sin"! W cos y sin (B cos y) dp 
e/7o 


_ B } sin”*! Y cos (B cos y) dy 


Ni eJo 
n+l i 
4 2? wT + ) 
a? ae Fn+i (B). (70) 


- * 
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Substituting these various integrations (67), (68), (69), and (70) 
in (66), we have 


12 . n he” n jn 
= — 4 — 1)? — 2 sin? B — 1)? ———_— 
P 4e | d | nin + 2 sin’ BY (— 1) n+ 1! 


n-1 n 
*) 9 Vv y 
n }” = v I (5) 





+ 2 cos BY (— 1)2 32 92. soy? _ T nt (B) 
B? ; 
n+l n . 
n }-” B 2 , Vir D (5 sa 
+ 2 sin BY (— | )2 9242. a “ —~C~C~*# Jn (B)- 
B? 2 
(71) 
2rb . T 
n=2,4,6,---@ . B= » Is between 0 and 5: 


This result may be expessed in a power series by expanding the 
Bessel’s Functions by equation (6), Art. 120 of Byerly’s Fourier’s 
Series, giving 








n-1 —- 
2 RB 
— ‘ + | ifs n+ 1) 
n a 
B* ~~ Be P 
+ vege (2H 1) n+ 3 2196 n+ 1) n+ 3 me 
= \ 2 / 2 oii » ) ~“ > 3 : 
(72) 
n+] 
: b 
Fe ] (B) = a B ; _ 
’ = a fs) 2 (243) 
ata. . 2 }| > | 
b Be 
al ae ee es SO a ae a 
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Note that 


bo 


-6- --n— 


4 
74) 
(3e5ee em —]’ (74) 








and 


, {nr 
S G * 1) 2:4-G---m 
r ( +3) 185 en + 1 








Putting these values in equation (71) we obtain 


PL er 1, sintBo 
p= CE S| a’ 3a@+D 
RB? B 1 


L | ees Tien 
vee. ices) * wes Dare 





Bs 1 i 
— 3189 @+)mt3)n+ts)* sarin ( 
B 1 
1 2 (nt+1) (n4+3) 
BS 1 
T 39 (n+ 1) a +3) (n +5) 





_pilj 1 
+ Ome les 


5 LS SERIO 
283 (n + 1) (n + 3) (n+ 5) (n + 7) \ |’ 


where 
n=z=24.6,.... (76) 


Equation (76) may be further improved for purposes of calculation 
by expanding the trigonometric functions in power series and col- 
lecting the terms. For this purpose 


sim B _ ~~ cos 2B _ 2?B4 24 B® 26 BS 


2 4 :- «2 © 
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-.F FF ” 
csB=1—5, +4, -@ + teens ' (78) 
4 6 
B sin B = BY te meee (79) 


Equations (77), (78) and (79) substituted in (76) will give 


Ir, n fn , 
pete An, Bn (B), (80) 


where F,, (B) is a polynomial in B®, B?, Bt, etc., where the coefficients 
of the several powers of B are contained in the table of page 225. 

In this table the bottom row of terms gives the coefficients of the 
powers of B, when the summation indicated in (80) is performed 
with n = 2,4,6----o. The various terms in the columns were 
employed in obtaining the last row by addition. 

The coefficient of B” is not contained in the table, because of its 
numerous terms, but its value when summed up is 


255n* + 6084n3 + 51396n? + 177264n + 193536 


10! (1+ 1) (n+ 3) (n+ 5) (n +7) (n+ 9) 
Substituting the values of the coefficients multiplied by the corre- 
sponding powers of B and summing up as indicated in equation (80), 
we obtain for the power the expression 





i ee We eiihoa | 
ie” 160 3780 | 56700 93555 ' 
_ 4) BY B Bo OB |] 
/ 1120 6480 83160 77395500 ' 
i) ae ae _. Se Ore ( —¢g 
Tt Hy 45360 ~ 24960960 * 6'34720 , fe 


This equation gives the average power radiated in the aérial hemisphere 
from the flat-top of the antenna regarded as a separate radiator with the 
distribution that it has under the fundamental assumptions of the problem. 
The current is to be measured in absolute electrostatic units, and the power 


is in ergs per second. 
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In this equation B 2b 


It remains to find how this power is modified by the mutual effect 
consisting of the interference between the waves emitted from the 
vertical portion of the antenna and the waves emitted from the hori- 
zontal part. This is the subject matter of Part ITI. 
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PART III. 


THe Murua. Term In PowER DETERMINATION. 


16. The Trigonometric Relations.— In Section 14, equation 
(56), it has been shown that the power radiated through an element 
of surface consists of three terms in the form 


dp = Pi (Ey Hy + Ey Hz + 2 cos a Ey Hy) dS. 


The first two of these terms we have already discussed. Putting 
in the values of EK, and Hy from equations (19) and (55) the remaining 
power term, which we have for convenience called mutual power, he- 
comes in the time average 


- I? cosadS . Az § 
dp = 


sin « ecosW sin B — sin (B cos y) 


me ro sin 6 sin ro | 


cos B cos (A cos #) — sin B cos8@ sin (A cos 8) — cos G . (82) 

In forming this equation we have multiplied the expression for Eg 
of eq. (19) by the expression for Hy, eq. (55). The product so obtained 
contains terms involving sin 7 cos r plus terms involving cos’r. The 
time average of the sin tr cos7 terms is zero; while the time average 
of cos? z is $; these facts have been used in forming (82). 

To be able to integrate equation (82) we must replace a, z, y and 
dS by their values in terms of 6, @ and 7%. By Fig. 3, 


= Tg COS 6, (83) 


a 
| 


dS = ro sin 6 d6 do. (84) 


In the spherical triangle of Figure 10, a is represented, as defined, as 
the angle between 6 and y, while opposite to a the side is 7/2. The 
important trigonometric relation in a spherical triangle is as follows: 

I. The cosine of any side is equal to the product of the cosines of 
the two other sides plus the continued product of the sines of these 
sides and the cosine of the included angle. 
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By this proposition, referring to Figure 10, we see that 


FIGURE 10. 


By the same proposition 
T o . 
cos 5 = cos 6 cos Y + sin @ sin y cos a; 


cos 6 cos y 
08 Oo ™ sin @ sin?” (86) 
or 
cosa —s- cos 8 cos (87) 


sin sin 6 sin? yp’ 
and by (85) this becomes 


cos a cos 8 cos d - 
—— wr mr (SS) 
sin y 1 — sin*@ cos* @ 


17. Integration for Mutual Power. — Now substituting the 
trigonometric relations (83), (84), (85), (88) into equation (82), we 
obtain the following integral expression for the time average of the 
mutual power radiated through the aérial hemisphere: 











T — 2 
cos Y = cos 5 cos 8 + sin 5 sin @ cos @ 
= sin @cos @. (85) 
6 
a 
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2 °7/2 ; 
p= | d6 sin (A cos 8) 


CT 


\ 


sin B cos 6 sin (A cos 6) — cos G 


cos B cos (A cos 9) — 


1 — sin’ @ cos? @ 


— cos @ sin 6 sin B is sana I (S89) 


» | — sin? 6 cos? @ 


{° cos ¢ sin (B sin @ cos ¢) do 
cos 6 olan cus 


This is a very complicated expression involving the integral of an 
integral. 
We shall first proceed to perform the integration with respect to ¢. 





Let V = | cos @ sin (B sin 8 cos ¢) de (90) 


1 — sin’ @ cos? ¢ 


and break the integral into the sum of two integrals thus: 


y= [+ i 


By a change of variable in the second of these two integrals by 
replacing ¢ by ¢’ + 7, we find that the integrand is unchanged, while 
the limits become 0 and 7, so we may write 


— [" cos ¢ sin (B sin 8 cos o) dp (91) 


1 — sin? 6 cos? @ 


Again decomposing this into the sum of two integrals we have 


y =2 |} oe [7 | (92) 


and changing the variable in the second integral by putting 
o = t — ©’, the second integral becomes 


[" 7 [° — dq’ (— cos o’) ‘a sin (B sin 6 cos ¢’) 
eJx/2 ‘ x /2 


1 — sin? 6 cos? ¢’ 





b 








230 PIERCE. 


which by dropping the primes and substituting in (92) and (91) gives 





va4 "sont sin (B sin 6 cos @) do 








93 
1 — sin? 6 cos’ (93) 
Now expanding in series as follows: 
; ° . B3 sin? @ 3 
sin (B sin @ cos ¢) = B sin 6 cos ¢ — a ? rN 
B’sin?6 co’ pd 
5! ' 


and 
1 


—s>——— = 1 + sin’ 6 cos? @ + sin*6@ cost‘é +..... (93a) 
1 — sin? 0 cos? 


and by taking the product of these two series we obtain 
°*x/2 
v= | di| B sin 800s? 
De Cs 4 
+ ) B 31 ¢ sin 6 cost d 
B? 


be 
—— p | cin? @coté 
s } Dp. ( 


+) B- 
a ceeseee i; (94) 


Integrating (94) by formula 483 of B. O. Peirce’s Tables, we obtain 


V = 2 E B sin @ 





1-33 \ B 7 
T 34 } Bs, ( sin” 6 
1-3-0 | i, 
. 2.4-6 ) Bs, T Fi ‘ho Y 
1°-3°35°7 \ Be B> Bi } ae 
v 2-4-6°S | b ~ 3} Bt om val ( sin! 6 


ee (95) 


We shall next proceed to perform the second integration with 
respect to @ indicated in (89). For abbreviation let us write 
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wo [~ tee 4 ( ees 


1 — sin? 6 cos? @ 1 — sin? 0 cos? 








by reasoning similar to the above. Expanding the denominator by 
(93a), we have 


W=4 ba “dé cos? 1 + sin? 6 cos? ¢ + sin‘ 6 cos* @ + - at 


a 


Cae ore ss : fi 
4+ 5. sin e+} 5. - > sin §+- ...-.-- ( (96) 





= 27 


oer ° 


nN) — 


(If we need it, this integral can be obtained by direct integration in 
the form 
\ l ( 


P is Mie ; 
W = on ( cos@(1-+ cos @) J 





but the expanded form is more useful for our purpose). 
Now substituting (95) and (96) in (89) we obtain 





*7/2 ' 
} d@ sin 6 sin (A cos 8) } cos B cos (A cos 6) 
— sin B cos @sin (A cos 8) — cosG 


l ' : 
: (B — sin B) sin 6 


3 . 
rae. B- ; — sin B) sin* sin® 


| 





- 3 5 
_ts 5 (p_ + % —sin B) sin’ 6 


2-4-6 3! 


davsuncews | (97) 


To evaluate this expression we must obtain the following integrals: 


T 
tie sin (2A cos @) | 
i= | dé sin” 6 = a : (98) 


°r/2 
[, = | d@ sin" 6 cos 6 sin? (A cos 8), (99) 
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*7/2 
I, = | dé sin® 6 sin (A cos @), (100) 


where n = 2, 4, 6, 8,....... 

I; is the simplest of these integrals and will be considered first. By 
expanding sin (A cos @) in series we have 
{~ ; A’ cos? @ , A® cos’ 6 


- 
: dé sin 9) A cos? — ——.,— — mm sees 


I, = 


which by Byerly Int. Calc., Art. 99, Ex. 2, may be integrated in 
Gamma Functions as follows: 











_ AB 2 
oe ———— ~*~ = 
ee ee, ar (*+3 41) 
2 
gst @)0 ns) 
T sy 
oo. 





Werrerrre eT (101) 


If we note that 


_(n+3 n+ 3 n+ 1 n+ 1 
r(y dF 2 2 ey) 





we obtain 


f {3 2 45 929! 


S@ O41 3! @+D@+3)' S@tD@ta@etrs 
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_ A } _ a. Ap A’ anes at 
n+l1 / 3(n+3) 35:°3(n+3) (n+ 5) 
AS he nas i. (102) 
7°53 (n+ 3) (n+ 5) (n+ 7) \ 7 
In like manner 
| 2A (24)? (2.4)! 
| h@® 37 Th: eo ae oe an 
: Z2(n+ 1) / 3 (n+ 3) 3 (r+ 3) (n+ 5) 
4 . 
: (2.4)° | cae 
| vee © (103) 
Z:a°3 neha Dash \ — 
: Now taking up integral I, from equation (99), let us write it 
| rf 1 — cos (24 cos 6) | 
: I, = | do sin" 6 cos 6 ) = = Rg 
J o / 2 \ 
and expanding cos (2 A cos @) in series, obtain 
] = \ (24)? cos? 6 = (2a)? cos* 6 } 
i= > al sin 9 cos 6 » 21 — 4 ie ilidadhdidid (- 


This equation, integrated in Gamma Functions between the limits 


rfat+l., 2)) 
= ( ; * 


0 and 7/2 gives 








l 
h= 3/5 | 
Z| 2. ar("y"+1) 
q 2 
n+1\_ ; 
watt 2 re 
ar geome “>< teks Sn : 
ar ( 9 +1) i 
_ 24 ),_247 __, 248 em, 
~ (n+1)(n+3) ) 3(n +5) 35-3(n+ 5) (n +7) y 


(104) 


Employing the values of i, Jz, Zz found in equations (103), (104) 
and (102) we may write the expression for the mutual power in the 


integrated form 
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, tila 
& 
| 
a. 
= 
SS 
—— 
IX 

= |b 

or to 

, oe 
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5+: 
1-3 BR 6. C\ AL, (GAP (24)* 
+54 (B- 31 sin B) P c T 9 


1°3°5 BS B A + 
+59 (F- sit 5: 7 8inB) [1 3-9 +5 


3°7 5.: je 





~ 3 5) 9 A2 ), 4 
L125 (p_ BP sna) Ly 24" 2a 
2-4-6 3! 5! 


3:11 a°3° 
(2.A)° 


11:13 


7-5-3-11-13-15 7 | 


a 
Af 
9+ 53. 


46 
7°9°3°9- 11-15 


eel 


9-11 


yd 


(104) 
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If now we recall that G = A + B, it will be seen that the equation 
(104) is entirely in terms of A and B and J. 

For purpose of computation it is found advisable to expand all of 
the trigonometrical expressions in power series and then perform 
with them the indicated operations. This was done with considerable 
labor and gave the following expression for mutual power: 


p= = Ec; 0166 Bt — .00404 B® + .000390 BS — 
0000144 BY +... 

+ 43 0083 B? — .00480 B® + .000729 B? — 
0000486 B® +: - 


. 
CV 


+a — .00433 Bt + .00104 B® — .000102 BS + 


0000051 BY —--. 


. r= 


-+- A® — .00127 B® + .000741 B° — .000106 B? + 


0000073 B? —--- : 

+ 4°} 000404 B‘ — .000101 B® + .0000101 BS — 
0000005 Bi + 
Se } (105) 


This equation gives the time average of the power radiated in the aérial 
hemisphere by the mutual effect of the fields from both parts of the antenna 
and is the correction to be added to the power radiated by the two parts, 
estimated as independent of one another. The current T is in absolute 
c.g.8. electrostatic units, and the power is in ergs per second. 


18. Summation of Flat-top Power and Mutual Power. 
We have obtained in equation (81) the time average of flat-top radi- 
ated power, and in equation (105) the time average of mutual radiated 
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power. If we replace the k of (81) by its value in terms of A, the two 
expressions may be added together. At the time of the addition 
we shall reduce the units to the practical system by multiplying the 
right hand sides of both power equations by 30 times the velocity of 
light in centimeters per second (i. e. by 30 c), and obtain 


_ 
5 = 60 P | 05 Bt — .00985 B® + .000849 Bs — 


0000356 BY... 
+ 43 0083 B® — .00480 B® + .000729 B7 — 
0000486 B® +. 
— 01148 B* — .00227 B® + .000198 B® — 
00000953 BY +... 
~ 45 ) 00127 B® — .000741 B®’ + .000106 B7 — 
0000073 B® +... 
+ 4° 00111 Bt — .00014 B® + .000019 BS—...... 
rae | (106) 


oo 


This 1s the total power contribution of the flat top by virtue of its indi- 
vidual and mutual action. The power is in watts, and the current I is 


in amperes. 
Certain Tables computed in the next Part of this communication 
make calculations with this series comparatively simple. 
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PART IV. 
COMPUTATIONS OF RADIATION RESISTANCE. 
19. Equation for Radiation Resistance. — If 


a = length of vertical part in meters, 

b = length of horizontal part in meters, 

\o= the natural wavelength of the antenna in meters, 

\ = the wavelength in meters of the antenna as loaded with inductance 
at its base, 


Pra 
A= : 
r 
Pah 
Bm —. 
To 
q = \? 


we may obtain the radiation resistance of the antenna by dividing the 
power radiated by the mean square of the current at the base of the 
antenna. This mean square current at the base of the antenna Is by 
(5) 
_ F sin’ (g/2) 
. 2 

Performing this division as to the flat-top power employing equation 
(106) and adding the result to the radiation resistance for the vertical 
portion as given in equation (44) we obtain for the total radiation 
resistance of the antenna the equation 


l \ 


aint (9/2) ) R, — Re cos = Rs sin q + 


Ro 2 


"> A? -. rsA® — ry AS — 7,5 + rel +-.-- ‘ (107) 
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This is Radiation Resistance in Ohms, where 





— bet we At? oy Ot 4e_.... | 
ittiees ) 3!2 war 5!4 i ie 7'6 ea) \ 
(2 +2-4., £42'-6,,,,, 
R, = 15 ) 319 (2.4) — 5'4 (21) + 
6 toa. (ZA) —-:-. ! 
4.0 \ 
ge NS PTR KI yg FHT —TF ya 
Re = 15) Gig 4) gp PAP 
oe ) 
re 
87 
ry = 120 05 B* — .00985 B® + .000849 BS — .0000356 B® +---- 
r, = 120 0083 B? — .00480 B® + .000729 B? — 
0000486 B® +.--- 


ry = 120 ) 01148 B‘ — .00227 B* + .000198 B’ — 


00000953 BY +- 
rs = 120 00127 B® — .000741 B> + .000106 B? — 
0000073 B® + 
m= 190 00111 Bt — 00014 B® + .000019 BS—---- ) 


(108) 


20. Tables of Coefficients of Radiation Resistance. — 
There follow in Tables I and II the values of the coefficients R;, Re, Rs, 
Te, 3, T4, 15, Tg for various values of A and B respectively. These tables 
have been computed by the equations (108). 
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TABLE I. 


239 . 





























COEFFICIENTS FR, Re, AND R3. 
2A \/4a R, Re R; 
. 31.416 04998 049919 002498 
3 15.70 19971 19870 01994 
3 10.47 44848 44344 06700 
4 7.85 79521 78107 1579 
5 6.28 1.2383 1.20634 3060 
6 5.236 1.7759 | 1.6969 5241 
PS 4.488 2 4055 | 2 2602 8232 
8 3.927 3.1240 | 2.8786 1.2137 
9 3.491 3.9290 | 3.5403 1.696 
1.0 3.141 4.8165 4.2315 2.300 
1.1 2 854 5.7837 | 4.9383 3.009 
1.2 ? 616 6. 8232 | 5.6442 3.823 
1.4 2.241 9.150 | 7.000 | 5.90 
1.5 ? O92 10.3392 | 7.611 6.999 
1.6 1.962 11.64 1 $8.15 | 8.35 
1.732 1.812 13.415 | 8.798 | 10.113 
1.8 1.743 14.40 | 9.10 11.20 
2.00 1.570 17.241 | 9.550 14.354 
2.20 1.427 20.15 | 9.55 17.80 
2 . 236 ~ 1.408 20.778 9 508 18.470 
2.40 1.307 3. | 9 00 21.42 
? .60 1.207 26.37 7.90 25.20 
2.642 1.189 27.053 | 7.60 | 25.927 
2.80 1.121 29 40 | 6.22 ' 29.05 
3.141 | 1.000 34.45 : 2.12 35.64 
TABLE II. 
(COEFFICIENTS fe, r3, ETC. 
| | 
B /4b re rs | r4 rs | rs 
1.4 | 1.112 15.8 686 | 3.55 055 | .481 
:.2 | 1.31 9.32 OSD | 2 ta 107 | 234 
| 
1.0 1.57 4.92 498 | 1.18 075 118 
8 1.96 2.16 340 | 513 051 050 
6 2.61 73 171 158 026 | .016 
| | 
4 | 3.93 15 057 085 009 | 004 
3 | 7 009 008 | 002 001 .0002 
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21. Curves of Resistance Due to Radiation from the Flat- 
top. — We shall now proceed to discuss the curves of radiation re- 
sistance of variously proportioned antennae when employed at various 
wavelengths relative to the natural wavelength. As preliminary, 
the resistance due to radiation from the flat-topped portion of the 
antennae is first computed. The equation for this is the summation 
of terms in (107) containing the small r’s as factors; that is, 


| ; / 
Ro — So ) ro A" + re Al? — ry? — rs? + reAl® + Ree - (109) 
“gin? (q/2) | \ 
due to 
filat-top 


in which 


2ra 
A= “% 
ah 
B 
q= ™ =2(A+B)- 


Since the coefficients (small r’s) are functions of B only, as given in 
Table II, it follows that when A and B are given, the value of the flat- 
top R may be computed. The results of the computations for various 
values of A and B are plotted in Figure 11. 

In this figure values of B are the abscissae, while the flat-top 
resistances in ohms are ordinates. The separate curves numbered 
1, .2, .3, ete. to .9 are for values of A = .1, .2, .3, ete. to .9. 

The outside end-points of these several curves, through which a 
limiting curve is drawn, are determined by the equality of the im- 
pressed wavelength and the natural wavelength of the antenna 9; 
that is, by the value of A + B = 7/2, which is the largest value 
A+B can have for the fundamental oscillation of the antenna. 


22. Curves of Total Radiation Resistance.— The next 
step consists in computing the radiation resistance of the vertical 
portion of the antenna, using the first three terms of equation (107), 
and employing a large number of values of A and B. To these values 
of resistance due to the vertical portion of the antenna the correspond- 
ing resistance of the flat-top are added so as to give the total resistance 


3 
; 








? em J 


Resistance of Top 
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of the antenna for various values of A and B. Curves of resistance 
for various values of A+ B are then plotted in Figure 12, with 
values of B as abscissae and values of resistance as ordinates. Figure 


1.6 


1.4 


LA 


_ 
© 


- 





Ps 4 6 8 im i 1.4 1.6 
B 


FIGURE 11. Resistance of horizontal top portion of antenna plotted against 
values of B. The separate curves numbered «1, *2, *3, ete. to «9 are for values 
of A = «1, «2, °3, ete. to «9. 


13 is an enlarged view of some of the curves that are on too small a 
scale to read in Figure 12. Then to make the family of curves more 
useful for ready reference a series of curves are drawn through all 
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' 
’ 
NY 
e 
oh x 6 Re) 1.0 1.2 1.4 
Figure 12. Total radiation resistance plotte] azainst values of B. The 
separate curves through the origin are for designated values of y. Separate 
curves not passing through origin are for different values of A + B. 
° 
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the points which have a common ratio of length of flat-top to length 
of total antenna. This ratio is designated by y, where 


B b 


77a ca (110) 


with b 


a = length of vertical part of antenna. 


length of flat-top 


These y-curves all pass through the origin. 
Next asa final step the curves of Figure 14 are taken from the curves 


RESISTANCE. 





4 
B 


Figure 13. Enlarged view of some of the curves of Figure 12. 


of Figures 12 and 13 with the new set of parameters. These curves of 
Figure 14 are the final curves of total radiation resistance, and are in 
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Figure 14. Total Radiation resistance plotted against \/Ao. The separate 
curves marked 0, *2, «3, etc. are for values of y = 0, ¢2, «3, ete. 
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terms of the ratio of the wavelength employed to the natural wave- 
length (that is \/Xo) and the ratio of the length of flat-top to total 
length of antenna (that isy). Figure 15 is merely a magnified view of 
certain of the curves that are too small to read on Figure 14. 


tf 


1.4 


1.2 


1.0 


RESIS TANCE 





2 4 6 8 10 12 14 16 
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FiGURE 15. Extension of curyes of Figure 14 to larger values of A/Xo. 


23. Total Radiation Resistance of a Straight Vertical 
Antenna at Various Wavelengths Obtained by Inductance at 
the Base.—.As an example, let it be required to find the total radiation 
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TABLE III. 


RESISTANCE OF A STRAIGHT VERTICAL ANTENNA FOR DIFFERENT VALUES OF 
WAVELENGTH OBTAINED BY INDUCTANCE AT THE BASE. 











d/o fe Radiation 
Ratio of Radiation Resistance in 
Wavelength to Resistance in Ohms 
on Computed by | computed on 
Present Theory 
1.00 a 36.97 98.7 
1.12 26.40 ie 
1.21 21.70 607.3 
1.31 | 17.65 D420 
1.45 14.28 48.2 
1.57 11.62 40.0 
1.74 9.10 32.6 
1.97 6.92 23.4 
2.24 9.19 19.4 
2.62 3.78 14.4 
3.14 2.98 10.0 
$95 1.65 6.40 
2.26 90 3.60 
7.85 30 1.16 
15.70 O82 40 
31.42 Ol .10 











resistance of a straight vertical antenna for various wavelengths 
obtained by adding various inductances at the base. For this case 
y = 0, and from the y = O curve of Figures 14 and 15 R may be 
directly read. The values which were used in plotting this curve are 
given in Table III, where they are compared with the corresponding 
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values computed on the assumption that the oscillator is a Hertzian 
doublet. This latter assumption gives 


ra 


{= 160 
R 16 x2 


It is seen that the departure of the present theory from the 
doublet theory is very large for the straight vertical antenna, as 
should be expected. 

It should be noted that the first value in the column of resistances 
computed by the present theory agrees with the value for this case 
computed by Abraham in the work cited in the introduction. This 
one value is the only value arrived at by Abraham for the fundamental 
oscillation, and is the case of a straight vertical antenna oscillating 
with its natural frequency. Abraham’s other computed values are 
for the harmonic vibrations with more than one loop of potential 
always without loading the antenna by inductance, and without any 
flat-top extension of the antenna. 


24. Comparison of Computations on the Present Theory 
with Dr. Austin’s Values for the Battleship ‘‘ Maine.’’—F igure 
16 gives the Radiation Resistance of the Antenna of the Battleship 
“Maine” as computed by the present Theory in comparison with 
Dr. Austin’s measured values of the total resistance of this antenna, 
and in comparison with values computed on the doublet theory of 
Hertz. The black dots of Figure 16 are Dr. Austin’s observed values. 
The heavy line was obtained by computation by the present theory, 
and the weaker line, by computation regarding the antenna as a 
doublet of half-length equal to the vertical height of the antenna. 

It is seen that the departure between the present theory and the 
doublet theory is not so great as in the case of the straight vertical 
antenna, for the reason that the doublet theory becomes more nearly 
correct as the half-length of the oscillator becomes small in comparison 
with the wavelength. 

Neither of the theories gives a rising value of the resistance with 
increase of wavelength, and, as Dr. Austin has pointed out, his rising 
values for long waves are probably not due to radiation from the 
antenna but possibly to dielectric hysteresis in the ground beneath 
the flat-top. 

I do not give more extended comparisons with experimental values 
at the present time, because I am now making some experiments to 
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see how much reliance may be placed in antenna resistance measure- 
ments made by buzzer methods of excitation in comparison with 
measurements made by excitation with gaseous oscillators and other 
methods of continuous excitation. 





400 800 1200 1600 2000 2200 
» 


as 

FicgurE 16. Total radiation resistance versus wavelength for the Antenna 

of the Battleship ‘‘Maine.”’ Black dots are Dr. Austin’s observed values; 

heavy line, computations by present theory; light line, computations by 
doublet theory. 


25. Example of Different Methods of Constructing an 
Antenna that Will Have a Specified Resistance for a Given 
Wavelength.— Let it be required to construct an antenna that will 
have a given resistance (4 ohms, say) for a given wavelength (2000 
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meters, say). To solve this problem, it is only necessary to look up 
the four ohm point on the different y-curves of Figures 14 or 15, and 
find the corresponding value of \/Xo. We can then find the Xo of the 
antenna, since \ is given. Dividing the \» by 4 we obtain the total 
length of antenna. The value of y gives the fractional part of this 
length which is to be horizontal. The complete result is tabulated 


in Table IV. 
TABLE IV. 


CONSTANTS OF THE DIFFERENT ANTENNAE THAT HAVE 4 Ounms RESISTANCE AT 
2000 METERS. 














| Total Vertical Horizontal Intensity 
, r/r \ | Length Length Length Factor 
} / NO NO in Horizontal 
Meters Meters Meters Plane 
.s | 1.075 | 1860 | 465 93 .0 372 2d 
| 
7 | 1.389 | 14385 | 359 107.7 251.3 300 
.6 1.67 1198 | 299 119.6 179.4 310 
| 
.o | 1.94 1030 | 258 129.0 129.0 312 
| 
4 | 2.18 916 | 229 137.4 91.6 313 
3 | 2.32 | 861 ) 215 150.5 64.5 314 
a 2.44 S20 | POD 164.0 41.0 10 
0 | 2.52 793 | 198 198.0 00.0 320 











The question as to which of these antenna to choose for the given 
purpose is now chiefly a problem in economics. The economic ques- 
tion is, which, for example, is cheaper: Two poles or towers 93 meters 
high and 372 meters apart, or one tower 198 meters high? This 
of course pre-supposes that it is designed to use a flat-top antenna 
instead of some other type, such as an umbrella. 

The problem is, however, not wholly economic because the lower 
antenna would be preferable as a receiving antenna on account of its 
weaker response to atmospheric disturbances. There is also the 
further question as to which of the tabulated antennae will give the 
greatest vertical intensity of electric and magnetic force on the horizon 
at a distant receiving station. This is the subject matter of the 
next Part (Part V). 
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PART V. 


FiELD INTENSITIES AND SUMMARY. : 
- " ’ 

26. The Electric and Magnetic Intensities at a Distant Point 
in the Horizontal Plane. — Equation. (19) wits the valnes‘of the 
electric and magnetic intensities at a distant point due to the vertical 
portion of the antenna. If we replace I of that equation by its value 
in terms of Ip from equation (6), and make cos 6 = 0, we have the 
intensities in the horizontal plane in terms of Io, which is the amplitude 
of the current at the base of the antenna. This gives 


2] . cos B — cos G| 
= 0 an ny 
Eo = Hg = — cos — (et — 1) To 
CTo r a 
2d 





(111) 





The quantities outside the square brackets are constant for a given 
distance 7) and a given amplitude of transmitting current Ip. The 
relative intensities are therefore determined by the factor in the square 
brackets, which we may designate by 


_ COs B- cos G 


Xx — 


— 
— 


2d 


(112) 


Using the values of B, G, given in equation (20) and the value of y in 
(110), this equation (112) becomes 


To Tg 
COS Y Dy = ooo >» 
, en A 


= 
sin 


Xo (113) 
2d 


This quantity X we shall call “The Intensity Factor in the Hori- 
zontal Plane.” It is to be noted that the electric and magnetic 


intensities in the horizon plane are not effected by radiation from the 
flat-top; for, by equation (55), the field intensities from the flat-top 
are zero for z = 0; that is, all over the horizontal plane through the 
origin. 
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Intensity Factor in Horizontal Plane 
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Pr, bad — . q . . . 
In Figure 17 the Intensity Factor in the Horizontal Plane is plotted 
lor various values of y and various values of \/Ap. Taking from these 
curves the values of the intensity factors corresponding to the values 


1 2 3 4 5 6 7 8 S) 
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Figure 17. Relative intensity of vertical component of Electric Force in a 
horizontal plane at a given distance from various antennae and for a given 
amplitude of transmitting current. 


of y and X/X of Table IV we obtain the results in the last column of 
Table IV. It is seen that the intensity factor is slightly smaller for 
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the larger values of the relative length of flat-top. This diminished 
value of the intensity factor should be compensated by the use of a 
slightly larger transmitting current. The required increase of current 
may be easily computed by equation (111). 


27. Summary. —This paper contains a mathematical theory of 
the flat-top antenna. The process employed consists in the integra- 
tion of the effects of an aggregate of doublets assumed to be distributed 
along the antenna so as to give a current distribution described by 
equation (1) and illustrated in Figure 2. The electric and magnetic 
field intensity due to each of the doublets is determined by the Maxwell 
and Hertz Theories for all distant points in space. These field in- 
tensities are summed up for all the doublets with strict allowance for 
the differences of phase due to different doublets; the summation 
gives the resultant field intensities. Then by Poynting’s theorem 
the power radiated from the antenna through a distant hemisphere 
(bounded by the earth’s surface assumed plane) is computed by the 
integration of a number of intricate expressions. From the radiated 
power the radiation resistance is obtained by dividing by the mean 
square of the current at the base of the antenna. Tables of coeffi- 
cients for computing radiation resistance are given, and curves are 
plotted of the calculated values of radiation resistance for different 
ratios of the length of the flat-top to the total length of the antenna 
and for different relative wavelengths obtained by loading the antenna 
with inductance. Curves are also given for determining the relative 
electric and magnetic field intensities in the horizontal plane for 
differently proportioned antennae variously loaded. Various equa- 
tions developed in the treatment may find application to problems 
in the design of radiotelegraphic stations. Although this investiga- 
tion was undertaken in ignorance of a simple case investigated by 
Professor Max Abraham, by a similar fundamental method, his work 
was discovered early in the course of the treatment and served as a 
check on one of the resistance values here given. This paper may be 
regarded as an extension of the remarkable work of Professor Abraham. 
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